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I. INTRODUCTION 

The problem of meson dissociation in QCD matter is of great importance in relativistic 

heavy-ion collisions. Indeed one of the goals of such collisions is to create and understand the 

f] 

quark gluon plasma (QGP) [1], and for many of the possible QGP signatures, the dissociation 
of mesons in hot or cold matter is involved at some stage. Understanding the mechanisms 
of meson dissociation is crucial in order to establish a quantitative description of the QGP. 
It is therefore important to establish how a fast moving bound state is broken apart due to 
the presence of QCD matter. 

One of the important experimental signature of the QGP formation in heavy ion collisions 
is the suppressed production of high-p^ hadrons. In standard medium-induced energy loss 
calculations [2], it is assumed that the partons from the hard scattering hadronize outside 
of the medium, and therefore that only the energy loss of partons is responsible for jet 
quenching. However for heavy D and B mesons, whose formation time is less than that of 
pions, hadronization can happen in the medium, and the problem of meson dissociation in 
QCD matter becomes relevant. Perhaps this could help improve the description of heavy- 
meson suppression [3], which is underestimated in the standard picture [4], even after the 



fl. 



inclusion of collisional energy loss 

Another measurement which is not yet fully understood is the suppressed production of 
quarkonia Q]. Strong J/ty suppression has been observed in nucleus-nucleus collisions by 
a number of experiments, and various theoretical explanations have been proposed. The 
natural explanation invoked is the Debye screening |7|]: a bound state dissociates because 
the attractive force between its quark and antiquark is weakened by the color screening 
caused by the QGP constituents, thermalized quarks and gluons. 

However, J/\I/ suppression is also seen in hadron-nucleus collisions and it turns out that 
cold nuclear matter effects are crucial [8| (note that the mechanisms of heavy quarkonium 



production in the vacuum are already quite involved 



The suppression seen in nucleus- 



nucleus col 
effects 10, 



ision is likely due to the interplay between Debye screening, cold nuclear matter 



and other phenomena such as recombination [12| or another mechanism for 
meson dissociation: multiple scatterings in the QGP. Interestingly enough, it has recently 
been shown that for infinite-extend matter, multiple scatterings (collisional dissociation) are 



a more efficient way of dissociating bound states than Debye screening [13J] 
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Concerning cold nuclear matter effects, it is important to understand them in order 
to distinguish to what extend experimental observations in heavy-ion collisions are due to 
initial-state or final-state effects. In the case of quarkonium production, cold matter effects 
have their own interest, they are relevant for vector meson production in deep inelastic 
scattering off nuclei at high energies, where large gluon densities in the nucleus are probed 
[14J]. This is an important part of the physics program at a future electron-ion collider 15]. 

In this work we focus on the multiple scatterings. This is the natural mechanism for 
meson dissociation in cold nuclear matter, and the dominant one in hot matter for a long 
enough medium. Having in mind high-energy mesons, we will work within the eikonal 
approximation. This is suited to discuss forward meson production in hadronic collisions 
involving only cold matter, and high— px production in heavy-ion collisions creating hot 
matter. We describe the interaction of the mesons with the target nucleus or plasma with 
the McLerran-Venugopalan (MV) [16( and Gyulassy-Wang (GW) 17| models respectively. 
We find that multiple scatterings of mesons are very similar in cold and hot matter, they 
are controlled by the saturation momentum, of either the nucleus or the plasma. This scale 
determines whether the meson sees a dense or dilute gluon density. We investigate in details 
the differences between elastic and inelastic scatterings of the target. 

The case where the target scatters elastically has been well studied, it only involves the 
calculation of a two-point function. The situation where one allows the target to break up 
has not received much attention, and this mechanism is the focus of this paper. It involves 
the derivation of a four-point function, which is the main technical result of this paper. 

nnn 

While several four-point functions have been obtained in the literature [18|, H9|, 12Q[ , the one 
computed in this paper is new. A general algebraic derivation is explained, and in the large 
iV c -limit a diagrammatic derivation is also given. 

The plan of the paper is as follows. In Section II, we show how the survival probability 
of a meson in a plasma and the production of vector mesons in deep inelastic scattering are 
related to a two-dipole correlator, in the eikonal approximation. In Section III, we compute 
this correlator in the MV model for cold nuclear matter and in Section IV this is done in 
the GW model for hot QCD matter with an emphasis on the large N c limit. Section V 
is devoted to discussions on possible consequences for the dissociation of quarkonia in the 
QGP, the suppressed production of heavy-mesons in nucleus-nucleus collisions, and for the 
production of vector mesons in electron-ion collisions. Section VI concludes. 
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II. DIPOLE SCATTERING IN THE EIKONAL APPROXIMATION 

In this section we derive the scattering matrix element between two qq color singlet dipole 
states. This is relevant for the calculation of the survival probability of a meson in a plasma, 
or the production of vector mesons in deep inelastic scattering. In both cases, these processes 
involve the two wave functions describing the fluctuation of the initial and final states into 
qq dipoles, and the S matrix element between the dipole states. We review the formalism 
in the eikonal approximation. 

We shall use light-cone coordinates with the incoming particle being a right mover. Using 
light-cone perturbation theory (for an introduction see 2l|) and neglecting higher Fock 
components, we write the qq dipole wave function of a meson or virtual photon with tri- 
momentum P = (P + , P±) and polarization A : 

/d^p &cc ~ 

J^3^hh(P±^ z )—^ \p±,z,h,v,P±-p±,l - z,h,c) , (1) 

where p = (zP + ,pj_) denotes the momentum of the quark, z is the fraction of longitudinal 
momentum P + carried by the quark, c and c are color indices, and h and h are polarization 
indices. 

The normalization of the qq state is 

<P',A'|P,A> = (2ir?2P + 5®(P>-P)5^J2 / %?l^>*)| 2 > ( 2 ) 

hh 71 

where the photon wave function is calculated in QED and is not normalized to unity, while 
the wave function for the meson is calculated in various models and required to be normalized 
to unity, meaning: 

E/%?lA(Px,z)| 2 = l. (3) 

hh J 

The S— matrix element (P'| S \ P) is 

Sfi = W / (2tt) 3 (lo 3 ^^' Z '^ P± ' z ^ S dc(P P) S cd( p -P P'-P) , ( 4 ) 

where S q and S q correspond to the scattering of the quark and antiquark respectively, and d 
and c are color indices. In formula (j3J), the polarization indices of the wave function are kept 
implicit, spins are conserved during the eikonal interaction and are not relevant. In fact the 
only quantum number changing is the color. The transverse momenta of the partons also 
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changes however in the high-energy limit, when the partons propagate through the hot or 
cold matter, they have frozen transverse coordinates and the matrix element depends only 
of those transverse coordinates. Therefore it is convenient to introduce the wave function 
and scattering matrix element in a mixed representation z) with x± = x q ± — x q ±: 

0(Rl, z) = I d 2 x ± exp <^ -ix ± ■ ( p± ^ — P± ) \ tp(z, x±) , (5) 

J I V m q + m q J J 

SI(P - P) = 2n5(zP + - z'P + ) J d 2 x q± e^- p '^W dc (x g ±) , (6) 

where (p± — m ' m ^ m _ P± j is the transverse momentum conjugate to x±. The Fourier trans- 
formation ([5]) is defined such that the P + dependence in f only enters through z = p + /P + , 
and such that there is no residual P± dependence in <p(z, x±). In formula ([6]), the scattering 
of the quark is described by the fundamental Wilson line (see for instance 221 ]) 

W[A](x ± ) = Vexp (ig s J dz + T c A~(z + ,x ± U (7) 

where V denotes an ordering in z + . This matrix in color space is a function of the classical 
color field A~ (we work in the gauge A + = 0) describing the QCD matter. The properties 
of this color field will be discussed in more detail in the following section. The scattering of 
the antiquark is described by the Wilson line W^ d (xq±). 
This allows us to write 

S fl = 2tt 2P + 5(P' + - P + )M (8) 



M 



J ^V^xe^^y^^^Xi)^^^)^!,^) • (9) 



The coordinate variables x± and X± are the dipole size and center of mass respectively, they 
are defined in terms of the quark and antiquark coordinates x q ± and x q ± in the following 
way: 

m q x q± + m q Xg ± 

x± = x q± - x q± , X ± = —^-i — ■ H -^— . (10) 

m q + m q 

The dipole scattering matrix S qq is a trace of Wilson lines: 

S qq (x q± , x q± ) = i-Tr [W (x q± ) (x q± )] . (11) 

Introducing the transverse and longitudinal overlap functions between the qq color singlet 
states (for completeness we have reintroduced the spin indices) 

^(^) = \ E / £ E <M ^X^> ^) . ( 12 ) 

A=±l J hh 
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(a) 



(b) 



FIG. 1: Two relevant cases which fix the center of mass of the dipoles X± : quarkonia (a) and 
heavy mesons (b). The dipoles have the same center of mass in the amplitude and the conjugate 
amplitude because we are integrating over the meson transverse momentum in the final state (|14p . 

and 

^E^/(^Vk^)> (is) 

hh 

the survival probability of the meson is given by 

PT ' L = ~A j J^^ M ^ = l d2x ^ d2x '±^l L ( X ^ T fl L *( X '±) ( S qq( X q±, X q±) S qq( x/ q±, x/ q ±)) : 

(14) 

where A = j d 2 Xj_ is the cross-sectional area and (• • •) represents the medium average 
which we shall discuss in the following sections. x' q± and the quark and antiquark 

transverse coordinates in the conjugate amplitude: 

x' q± = X ± + -^—x' ± , x' q± = X ± - — 5— ^ . (15) 
H m q + mq H m q + mq 

We have assumed that <Tr [W F (x a± ) W F] {x b± )] Tr [W F (x' b± ) W F] {x' a± )] ) is independent 
of X± which we justify in the following calculations. 

The diffractive production of vector mesons in deep inelastic scattering also involves this 
dipole-dipole correlator, the cross section reads 

<7T,L = j d 2 X±d 2 x' ± $y^(x±)$y^*(x' ± ) j (fX ± ([1 - S q q(Xq ± ,Xq ± )}[l - S q q(x'q ± ,x' q± )]) , 

(16) 

with the appropriate overlap functions $y 7 between the photon and vector meson 
wave functions. Note that if we replace the correlator \S qg -(x q ±, Xq±)S q q( y x' q± , x' ql )} , by 
(Sqq(x q ±,Xq±)) ( y Sqq(x' q± ,x' q± )') , then we are only including scattering processes in which 
the target scatters elastically, this corresponds to exclusive production. However formula 
( |T6l) includes the possibility of the target dissociating. 
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The purpose of this paper is to calculate (S qq (x q ±, x q ±)S qq (x' q± , x' q ±)) , and to compare 
the result with the product of two-point functions {S qq (x q ±,x q ±)) (S qq (x^ ± ,x' q± )") . We shall 
compute the medium average in the MV model for cold matter, and in the GW model for 
hot matter, the calculation is essentially the same. The derivation of this 4-point function is 
the main result of this paper. As pictured in Fig. [I], we will be interested in two particular 
cases: 

• (a) the quarkonium or vector meson case m q = m q meaning: 

• (b) the heavy meson case m q 3> m q , this puts two Wilson lines at the same coordinate, 
which greatly simplifies the calculation: 

X±_ = x q± = x' q± . (18) 

Note that we shall only consider quantities integrated over P' ± , the transverse momentum of 
the meson in the final state. This is why the two center of masses in the amplitude and the 
conjugate amplitude are the same (see formula (TH]) ). However, the P' ± dependence could 
also be obtained from the results derived in the following. 

III. MULTIPLE SCATTERING OF A COLOR SINGLET DIPOLE IN THE 
MCLERRAN-VENUGOPALAN MODEL 

In the Color Glass Condensate framework, the low energy partons of a nuclear wave 
function, those relevant in high-energy processes, are described by classical color fields. The 
MV model Ji| is a model for the distribution of color charges which generate the field. It is 
a Gaussian distribution whose variance is the transverse color charge density squared along 
the projectile's path /i 2 (^ + ). The only parameter is the saturation momentum Q s , with Q 2 S 
proportional to the integrated color density squared. 

A. Introduction to the MV model 

The nuclear average of a function of Wilson lines / [A] reads 

(f[A)) = j Dp exp (- J d 2 xd^dz + Pc{Z+ 2 ^ Z ^ V) ^ f[A] , (19) 
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where the color charge p c and the field A~ obey the Yang-Mills equation 

-V 2 A-(z + ,x)=g sPc (z + ,x) . (20) 

The MV distribution is a Gaussian distribution, therefore one can compute any average by 
expanding the Wilson lines in powers of gsPc and using Wick's theorem. All correlators of 
p's can be written in terms of 

(p c (z + ,x)p d (z ,+ ,y)) = 5 cd 6(z + - z' + )5 {2 \x - y)p 2 {z + ) . (21) 

Note that we dropped the _L indices denoting transverse vectors, in order to get lighter 
expressions in the following. 

Inverting equation (1201 gives A~ in terms of p c : 

/f Ml, ik-x 

d 2 zG(x-z)p c (z\z) , G(x)= J (22) 

\k\>A QOD 

where G is the two-dimensional massless propagator. After expanding the Wilson lines and 
applying Wick's theorem, every contribution is a product of correlators like 



g 2 s (A-(x + ,x)A d (y + ,y)} = S cd 5(x + - y + )^ 2 (x + )gl j d 2 z G(x - z)G(y - z) (23) 

= 6 cd 5(x + -y + )p 2 (x + )L xy (24) 

times a trace of color matrices. The color algebra is the difficult part to deal with. 



B. The dipole-dipole correlator 

We now compute the following average 

(S q q(x,y)S q q(u,v)) . (25) 

Let's represent each W by a line along the z + direction, at a given transverse coordinate. 
Due to their respective z + ordering, W*8 are oriented to the right and W^s to the left. Due 
to the color structure, the lines are connected as shown in the left diagram of Fig. [2J Let's 
expand the Wilson lines in and use Wick's theorem. Graphically, every (AA) correlator 
(given by (T53J)) can be represented by a gluon link between two Wilson lines, at the relevant 
time z + , and the factor associated with it is p 2 (z + )L xy with x and y the transverse positions 
of the Wilson lines. This comes with a minus sign if the line is between two (anti) quarks. 
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FIG. 2: On the left is a representation of the correlator (I25p with an horizontal line for each Wilson 
line. Each vertical link corresponds to a (AA) correlator (|24p in the gsA expansion. The figure 
shows how, when using the Fierz identify (127j) . the color structure can be obtained in terms of the 
two coefficients a n and b n . Including diagrams in which a Wilson line is connected to itself, one 
gets formula ([28|) . 



A first class of diagrams, easy to deal with (and shown later in Fig.[6](a)), are those where 
a Wilson line is connected to itself. Then the contribution is (for instance) —Cf^ 2 L xx /2, 
it is color singlet and factorizes. The 1/2 is due to the z + ordering in a single Wilson line. 
Summing contributions with an arbitrary number n of such links (and diving by n\ in order 
not to overcount diagrams) yields 

^ g 2~ ^ (LxX-\-Lyy-\-L UU -\-L VV ) (26) 

The second class of diagrams are those where two different Wilson lines are connected, as 
shown in the figure. Let's simplify their color structure by using 

T t a j T ki = -j$u$ik ~ 2N$ij$ki (27) 

which is represented by the lower diagram in Fig. [2j In doing so, only two topologies can 
be obtained, the (x, y) and (u,v) color connections can be flipped into (x,v) and (u,y), or 
not, and any diagram is the sum of two contributions. 
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Therefore one can write 



(S q ,(x,y)S gq (u,v)) = ^f2 [ [N 2 c a n (z+,...,z+) + N c b n (4,---,4)] ■ (28) 

C „_n J 



n=0 



The diagrams have been classified by the number n of links they contain. At each order, 
the sum of diagrams is made of two contributions. The ones where the color structure is 
not changed (whose sum is denoted a n ) are multipled by N 2 and the ones where the color 
structure is flipped (whose sum is denoted b n ) are multipled by N c . This is illustrated in 



a n and bn can be obtained by iteration. The derivation is given in 







' a n -i \ 


U ) 




\ bn-l ) 



Fig. [2J The values o: 

the appendices A of jj] and 23j]. The sum of diagrams at order n is given by the sum of 
diagram at order n — 1, times the factor for the n-th link. There are six possibilities to add 
the n-th link on each of the two classes of diagrams. This gives 



(29) 



where the 2x2 matrix is given by 

u | {L xy + L UV )C F + ^- c F(x,y;u,v) -±F(x,v;u,y) ^ 

-~F{x, y; u, v) {L xv + L uy )C F + ^F(x, v; u, y) 

with 

F(x,y;u,v) = L xu -L xv +L yv -L yu = g A s J d 2 z[G{x-z)-G{y-z)][G{u-z)-G{v-z)} . (31) 

The problem has been reduced to finding the eigenvalue \± and eigenvectors of M. Indeed 
one has 

aH-(:)-(sH(«:«) - 

with a± and b± obtained from the eigenvectors of M. The results read: 

A± = (^j- ~ (L xy + L uv ) + ^(L xv + L uy ) + V\ «, v) ± ^v 7 ^ (33) 

\/A±F(x,u;y,v) F(x,y;u,v) 
a± = = , o± = =F ■= — (34) 

2Va n c Va 

4 

A = F 2 (x, u; y, v) + j^F{x, y; u, v)F(x, v; u, y) . (35) 
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Finally, resuming the contributions yields 



T 

(S q q(x,y)S q q(u,v)) = — 



N 2 \a+e lj2x+ + a_e^ A ' ) + N c (b + e^ x+ + 6_e M 



with 



/i 2 = J dz + iS(z + ) 
One can write the final result in the following form 



(36) 



(37) 



(S q q(x,y)Sgq(u,v)) = e 2 



^f[F(x-y)+F(u-v) 



= (Sqq(x,y))(S q q{u,v)) 



F(x,u;y,v) + VA _ F(x,y;u,v) \ 
2\/A N?\/A 



F(x, u; y,v)-VA F(x, y; u, v) \ 



2y/A N*VA 
It is given in terms of a single function 



e -^^F(x,u;y,v)+^^F(x,y;u,v) 



(38) 



F(x -y)= /x 2 (L xx + L yy - 2L xy ) = g^ 2 J d 2 z [G(x - z) - G(y - z)f 



Indeed also one has 



2/z 2 F(x, y\ u, v) = F(x — u) + F(y — v) — F(x — v) — F(y — u) . 



(39) 



(40) 



Note that for x — y or u = v, F(x, y; u, v ) = and one recovers the single dipole average. 

In the function F(r), the infrared cutoff Aqcd only enters through a logarithm as ex- 
pected. In the \r\AQCD^^ limit, one has 

1-Mk\r\) 



-F(r) 



2tt 
ff_ 9$C F 

4 47T 



dz + /i 2 (z + ) ) / (Ik 

Aqcd 



k 3 



dz + ^(z + ) lo. 



r 2 A2 
' IV QCD 



(41) 



This is the standard definition of the saturation scale in the MV model. In the following, 
we will neglect the logarithmic dependence of Q s , meaning 



F{r) 



\i F(x,y;u,v) 



-(x-y) ■ (u-v). 



(42) 



2C F v ' 2C F 

Finally, note that this model is easily generalized to hot matter. It is showed in Ref. [13j that 
by only keeping terms up to first order in k±/k z , the gluon propagator in light-cone gauge 
for a color charge propagating along the light cone in hot matter is 
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where y? D = \g 2 T 2 (N c + -f) and rf = ^=(1,0,0, -1). To generalize the MV model to hot 
matter, we treat the valence parton distribution as recoilless color sources which are localized 
along the light cone, and in the eikonal approximation, the only change of this model is to 
redefine G{x) as 

^-JtTT > ( 44 ) 



r) 2 P + /i^ ' 

and the infrared cutoff Aqcd is replaced by the Debye mass of the thermal plasma. In 
the next section, we will see our generalized MV model is equivalent to the GW model in 
light-cone coordinates in light-cone gauge. 



C. The quarkonium case 



Let us consider the situation where the two dipoles have identical center of mass: 



x=X+- y = X u=X V =X+— 

2 y 2 2 2 



(45) 



Then one finds 



F(x,y;u,v) 



2C, 



-r ■ r' F(x, u; y, v) 



Qs ( , /\2 

W r+r) 



(46) 



which yields 



(S qq -(r)S qq -(r>)) = e—i 



-%(r 2 +r'2) 



2VA' 



_( r I r 'Y /4+vA' r ■ r' \ 9i^AL 
— 2: — ■ — >—L — 1 _ |e 4 



JV C VA' 



+ i 



r + rT/4+vA' r ■ r' 



2VA 7 



\ e 4 1-1/JV2 



<£$_ T-T 

16 1-1/N'i 2 jv2_! 



(47) 



with 



16 7V C 2 



The large— N c limit gives: 



(^-(r)^-(r')) = e-- 



Qs r„2 , „/2 



(r 2 +r' 2 ) 



16 (rV) 2 / Q 2 (r + r^ 2 
1 ~ iV^r + r') 4 ; + 



I 16 ( r - r T _-gj(r-rT 

A^2 (r + r') 4 



(48) 



(49) 
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D. The heavy meson case 



We now consider the situation where both antiquarks are at the same position: 



x = X + r y = X u = X v = X + r' 



(50) 



Then one finds 

F(x,y;u,v) = F(x,u;y,v) 
which significantly simplifies the result: 



Q K-r' 



2C 



(S qq (r)S qq (r>)) = e 



-^(r 2 +r' 2 ) 



The large A^— limit is simply 



(51) 



(52) 



(S qq (r)S qq (r')) 



e 4 



Nf y + 2 rr ) + Nf 



(53) 
(54) 



IV. MULTIPLE SCATTERING OF A COLOR SINGLET DIPOLE IN THE 

GYUL ASSY- WANG MODEL 



In this section, we use the GW model 



171 ] to deal with the multiple scatterings of a 



color singlet dipole in hot QCD matter. The special case of heavy-meson dissociation due to 
multiple scattering is addressed in jsj , however the color structure is ignored in that analysis. 
In the following, we will show that the GW model gives the same results as MV model in 
terms of the saturation momentum Q s in the eikonal approximation. 



A. Introduction to the GW model 



In the GW model, the medium is modeled by an interaction Hamiltonian with iV — > oo 

scatterers: 



i=l ai=l 



X 



^(x)T a ^(x)^,(x) + ^\{x)T a *V t {x)^ q {x) = H q (t) + H q {t) . 



(55) 
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The screened potential 



Vi{x) = " e"" 01 *-* 1 , (56) 



or Vi(q) = e~ iq ' Zl in momentum space, is characterized by the Debye mass fio, and 

describes the medium in the situation /i£>A ^> 1, with A the mean free path of a single 
(anti) quark. Even though we are dealing with multiple scattering of a color singlet qq dipole 
in this paper, we still use the mean free path of a single (anti) quark and treat the quark and 
antiquark in the dipole wave function as two individual free particles. This approximation 
is good when the relevant length scale of the medium L < 77", the typical time scale in a 
meson [l3]. The assumption 1///D <C A means that the scatterers are independent of each 
other and, therefore, completely uncorrelated. 

The ensemble average over the transverse positions of the scatterers is defined as 

N „ , 2 

<•••>=£/ ^> ( 5 7) 

and after the medium average, each pair of color indices dj of the generators T 0,1 at the 
position Zi in the amplitude and/or the conjugate amplitude is identified as illustrated in 
Fig. [3(a) , which is equivalent to ( 1241) in the MV model. This equivalence enables us to 
calculate the color structure in the same algebraic way than showed in Sec. III. Instead, in 
this section we choose to calculate the dipole-dipole correlator in the large- iV c limit, which 

oo 

allows a diagrammatic analysis. We shall calculate the probability P = P^ in the 

n=0 

eikonal approximation with 

pW = i/§< |MW|2 > < (58) 

the contribution to P from those diagrams with n scattering centers scattering with the 
dipole in the amplitude and/or the conjugate amplitude. 

Among N scatterers, we have = u ^ ~ choices of the n scatterers which 

° ' JV n\(N— raj! n! 

scatter with the dipole in the amplitude and/or the conjugate amplitude. In the eikonal 
approximation, each choice of those n scatterers will give the same contribution to P^ n \ 
Therefore, let us assume that the dipole, moving along the + direction, is scattered by the 
first n scatterers which are located at z},i = 1, ■ • -,n respectively with zf > zf_ x . In the 
following, we enumerate the n scatterers by Zi and by Z\ > Zj we mean zf > z~j . Focusing 
on the color factor of each term contributing to P^ n \ there exists a Feynman diagram with 
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FIG. 3: The equivalence of the color factors between (a) and (b): after the medium average, each 
pair of color indices cjj of the generators T ai at the position Zi in the amplitude and/or the conjugate 
amplitude is identified. Therefore, it has the same color factor as the Feynman diagram (b). In 
the large N c limit, the amplitude of (b) is proportional to Cp in contrast with that in Fig. Hfb). 

the same color factor. An example is showed in Fig. [31 We will use this correspondence to 
evaluate P by a diagrammatic analysis at large N c . Up to 0(a 2 ), we will evaluate separately 
the following two cases: 

• single-scattering diagrams such as those shown in Fig. H](a). In this case, since the n 
scattering centers all scatter with the dipole by single scattering in the amplitude, they 
must also scatter with it in the conjugate amplitude. This corresponds to processes in 
which the target scatters inelastically. 

• double-scattering diagrams such as those shown in Fig. [3](a). In this case, if the 
scatterer Zi undergoes double scattering in the amplitude, there must be no scattering 
between it and the dipole in the conjugate amplitude and vice versa. This corresponds 
to processes in which the target scatters elastically. 

B. The evaluation of single-scattering diagrams 

The amplitude for the dipole to undergo n single scatterings with those n scatterers 
is a sum of 2 n corresponding diagrams since the gluon line from each scatterer can 
hook either on the quark or the antiquark line. If the gluon line from the scatterer Zi 
hooks on the quark or the antiquark line, the corresponding amplitude picks up a phase 
e l v qj - i ' x± = e l U-i-v x qj--*±) or e~ l ^r g±t ' Xj - = e* 9 - 1 -*' 12 ^ J-- 1 respectively, where q±i denotes the 
transverse momentum of the gluon line. In the mean time, the color matrix T Zi is put in the 
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C 



I) 



C 



3 



Zl Z2 Z3 Z4 Z4 Z3 Z2 Z1 

(a) 



Z1 Z2 Z 3 Z 4 
(b) 



FIG. 4: An example of single-scattering diagram of the leading order in N c : The contribution of 
(a) to P( n ) has the same color factor as the amplitude of (b), which is proportional to Cp in the 
large N c limit. 



quark part or antiquark part of the trace of the n color matrices accordingly. For example, 
if the scatterers z ai < ■ ■ • < z am and Zb l < • • ■ < z\ >n _ m scatter with the quark and antiquark 
respectively, we have Tr [T Zam ■ ■ ■ T Za ^T Zb i ■ ■ ■ T Zb "- m ] in the amplitude of the corresponding 
diagram. In general, we have 

• n—l r j9 

arqu 



n 



Airae 



-iq±i-z ±i 



Q±i + f*r> 



i-^q±i-x± T ai x 



- e -i-£*±i-*±T? H Sac.., 

+1 didi+i <-j<-i+l 



(59) 



with c\ = di, d n+ i = c n+ i and M = m q + mq. 

In fact, it is easy to show that we only have 2 n ~ 2 different traces of color matrices in 
which correspond respectively to the 2 n 2 different hookings of the gluon lines from 
the n — 2 scatterers in between z\ and z n on the quark and antiquark lines. Generally, we 
write the trace in the following form 



r^.(n,m) ^ rj r ]rpz n rpz am _ _ _ rpz ai rp Zl rpz bl _ _ _ y a &„_ m _2 ] 



(60) 



and the corresponding four diagrams give to Ms^ a contribution 



n-l 



( 1 \n— m— 2 r 

M (n,m) = (^l ) ^(n.m) / ^ x± ^ x± ) / J j 



d 2 q±j 
(2vr) 2 



x 



n 

i=i 



e i=i i=i 



(61) 



x e M 



J,Ii ) (< 
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FIG. 5: An example for 1/N 2 suppressed diagram of single-scatterings to p( n ): in (a) the gluon line 
from the scatterer Z2 hooks on different fermion lines in the amplitude and conjugate amplitude. It 
has the same color factor as the Feynman diagram (b), which is non-planar and 1/N 2 suppressed 



relative to Fig. H|b) 



2M- 



In the large N c limit, to evaluate 



p{n) _ n 



i r d 2 p' 



;|mW| 2 > 



(62) 



A J (2tt) 2 

we only need to take into account those diagrams with each pair of the n — 2 scatterers in 
between z\ and z n hooking on the same fermion line in the amplitude and conjugate ampli- 
tude. For example, in the case n = 4, the contribution to from the diagram Fig. 0Ja) 
has the same color factor as Fig. SJb) oc Cp, while the color factor of the contribution from 
Fig. G3^a) correspond to a non-planar diagram Fig. Mb) and, therefore, are 1/N 2 suppressed 
relative to Fig. H(b) 24|. This means that we only need to count in the 2 n ~ 2 products of 



M^' m ^ and their conjugates, which have a unique color factor 



Tr [T Zn ■ ■ ■ T Z1 ] Tr [T 21 • • • T Zn ] = C F + (-1) 

Therefore, for arbitrary n > 2, we have in the large N c limit 

1 



N 2 - 1 

1 c 

2 n N n 



CI 



(63) 



pW=-P 2 (Z + 2 X ) 



n-2 



n 



(64) 



where 



P 2 = j^(Lpa) 2 J d 2 x ± d 2 x' ± ^ fl (x ± )^Ux'_ 
d 2 q± 



n 2 



— e 



. mq 



^ J d 2 x ± d 2 x' ± ^%(x ± )^tW±)\Qt(^ ■ ^ 



-i^fq ± -x ± _ i^q^xj 



(65) 



/ \2 
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(a) Type I DS (b) Type II DS 

FIG. 6: Two types of double-scatterings for each scatterer in the amplitude: after the medium 
average the color factor of (a) is trivially Cp while the color factor of (b) depends on its relative 
position with respect to those for single-scatterings. If it is not in between any two of single- 
scatterings, it also has a color factor Cp. 



and 



Z = [ d 2 / ^.(^-VJ + e ^ q±<x± - x > ±) _ 2 \ (66) 

We have introduced the number density in the transverse plane ^ = Lp with L the 
plasma length and p the number density of scatterers, the cross section of the (anti) quark 
undergoing a single scattering a = A7r ^° F ; and the average scattering number x — j — Lpa. 
Finally, 

kj.KVki.1 

is the plasma saturation momentum squared, which is a characteristic property of the QCD 
medium. It is different from the traditional saturation momentum introduced in the previ- 
ous section (which characterizes the small— x part of a hadronic wave function), but when 



writting (1671) in terms of the gluon density per unit of transverse area in the plasma [251 ]. 
one finds the same expression than when writting the saturation scale in terms of the gluon 
density per unit of transverse area in a hadron. And indeed, one can check that formula 
(jUJ) is recovered in the |x_l|/zd *C 1 limit if one replaces pn with Aq C d m fl67)l - In terms of 
the traditional jet quenching parameter q, one has Q 2 S = hatqL. Note also that in formula 
(|65|) . we have assumed that \q±_ • 1. This is justified because the typical values are 

\q±\ ~ T while \x±\ ~ as, the meson size. We always assume a situation where Tas 1, 
where multiple scatterings are important. 



19 



C. The evaluation of double-scattering diagrams 

In all the diagrams, we only need to evaluate two types of double-scatterings as showed 
in Fig. El Each type I double-scatterings and its conjugate give a contribution —2% to P^ n \ 
This is easily read out from the corresponding diagrams such as the type I double-scattering 
at z 2 in Fig.[3](a). Those diagrams with at least one type II double-scattering in between two 
single-scatterings are 1/N% suppressed. For example, the contribution to corresponding 
to Fig. [3]^a) oc C\ is 1/AT? suppressed compared with that in Fig. @Ja) in the large N c limit. 
Otherwise, each type II double-scattering and its conjugate give to P^ a contribution 

Y + 2 X = f d 2 q ^ ( L 2 P ^\. 2 (e^ + e~^<) , (68) 



with 



Y^-\Ql{xl + x'l). (69) 

If the dipole interacts with the medium only through double-scatterings, that is, the target 
is not excited, we get the probability for the so-called elastic process 

f °° 1 

J m ^ miXmn - (70) 

= J rf 2 x ± rf 2 x' ± $^(x ± )$^(x / ± )e-^' ( ^ + ^ ) . 

D. The survival probability from inelastic multiple scatterings 

Now, we are ready to evaluate P^ el up to order 1/N 2 , the probability for the transition 
of an initial state tfi to a final sate (p / undergoing at least n > 2 single-scatterings with the 
medium excited. Let us assume that there are mj type I double-scatterings and rxiji type II 
double-scatterings among total n scatterings. First, we have = nl ^L n y — choices of 
the n scatterers. Then out of those n scatterers, we have C™~ mi possible choices of (n — mi) 
scatterers undergoing (n — mi — m n ) single-scatterings and ran type II double-scatterings. 
Since we only need to count in the diagrams with no type II double-scattering in between 
any two of all the (n — mi — mn) single-scatterings, we only have (mn + 1) choices of the 
m n scatterers undergoing type II double-scattering. 
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Therefore, for each group (n, m/, m.//), we have 



P (n,mi,mn) = C" 1 (m n + l)P 2 (-2 X ) mi (Y + 2 X ) mn (Z + 2x) n ~ mi ~ mu ~ 2 
n! 



(71) 



m/!(n — m/)! 



which leads to 



n— 2 n—2—rrijj 

EV~^ p{n,m u m u ) 
/ j 

mjj=0 mj=0 



1 n-2 

-V(m + i)y m z 

n! ^ 



n— m— 2 



m=0 



and finally 



inel 



n=2 



(n) = ^2 



[e* + (y _ Z _ i) e Y] 



(72) 



(73) 



We can now discuss the two situations of Fig. CD (we recall that x± = x q ± — x q ± and 

rp rfJ rf ^ 

x ± ~ x q± x g±J- 

• (a) In the quarkonium case: m q ~ m q and X± ~ (x q ± + x q ±)/2 = (x' q± + x' q - ± )/2, this 
leads to Z = — %-(#± — ^jj 2 an d 



PLl = ^2 / rf^±rfV ± $),(x ± )$J*(x 

( 16(q: ± -^) 2 c _i Q 2 (a;| _ ;r , ) 2 _ 2(x ± ■ x' ± ) 2 
\ (x ± +x' ± ) 4 (x ± +x' ± ) 2 



(x± + x' ± y 



(74) 



(b) In the heavy-meson case: m q 3> m q , and X± ~ = this leads to Z 



— %(^± — ^l) 2 and 



= J d 2 x ± d 2 x' ± <S> x fl (x ± )$ x f *(x' 



(75) 



In both cases, we obtain the same survival probability 



pX _,_ p\ 

el ' inel 



rf 2 x ± dV ± $) i (a; ± )^(a; / ± )( > S w -(x ± )^-(xl)) 



(76) 



as in the previous section, the differences enter through the saturation scales only. The 
infrared scale is Aqqd m the cold QCD matter case, and for hot QCD matter. 
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V. PHENOMENOLOGICAL IMPLICATIONS 

It is very interesting to see how the two different approaches, for cold and hot matter, yield 
the same result when expressed in terms of the appropriate saturation scale. This fact allows 
us to make general comments about bound state dissociation and vector meson production 
leaving all the medium dependence in only one parameter. From the calculations shown in 
previous sections we can easily see that in the large N c limit, the main contribution to the 
dipole-dipole correlator comes from elastic processes. The contribution coming from inelastic 
processes is 1/N% suppressed, and is therefore 10 % smaller in general. However, the different 
dependences of the subleading contributions on the size and orientation of the dipoles can 
play an important role in determining when these contributions become comparable to the 
one from the elastic case. For instance, we notice that the elastic correlator (as shown in 
(170]) ) is insensitive to the relative orientation of the dipoles but the inelastic contribution is 
not. 

First, it is worth noticing that the inelastic contribution goes to zero as Q\ when Q s goes 
to zero both for the heavy meson and the quarkonium case. In this limit which corresponds 
to turning off the interaction, elastic dissociation goes to zero only as Q^, therefore in order to 
get a large contribution from inelastic processes Q s should be large enough. This observation 
has to be taken with caution since Q s is present in all of the exponentials. A too large Q s 
induces a big suppression in the dipole-dipole correlator (elastic or inelastic) except for very 
small dipole sizes, which don't give a big contribution when the wave functions are included 

(Eq. (ED, (USD). 

With the above considerations in mind, we can easily notice that the terms in the inelastic 
part of the correlator with the exponential factor e~^ L ( x - L+x± ^ in Eqs. ( 1741) and (IT5|) don't 
give a large contribution, since they only become comparable with the elastic part in a 
region heavily suppressed by the exponential factor and the wave functions. Following 
this observation we turn our attention to the terms with the exponential factors of the 
form e~~ (x± ~ x± > and e~~ (x± ~ x - L > . Unlike the other terms, these exponentials are not 
necessarily small for large dipoles and don't exhibit the property of color transparency, 
which is characteristic of the elastic terms. As long as the two dipoles are aligned and 
similar in size, the first term in both expressions for the inelastic part can overcome the 
1/iVg suppression and become comparable with the elastic part. 
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FIG. 7: Survival probabilities for J ftp mesons, using "light-cone gaussian" wave functions and the 
dipole-dipole correlator 

More precise comparisons of the relative sizes for different dipole configurations require 
a more intricate numerical analysis which is left for future work, especially the result of 
the interplay between the size dependences of the meson wave functions and those of the 
dipole-dipole correlator should be investigated numerically. Still, let us make more specific 
comments concerning different physical situations where we expect our results to be relevant. 



A. Super-penetration of quarkonia 

As an example of a full numerical calculation, we considered the J/ ip survival probability. 
The light-cone wave functions of Ref . 26j| were used to obtain the probability plotted in Fig. [7] 
as a function of Q 2 S . Current measurements in heavy- ion collisions at RHIC suggest a value 
of Q 2 S ~ 1 GeV 2 for a gold nucleus and Q 2 S = qL ~ 20 GeV 2 for a plasma with a size of a few 
fermis for central collisions. At the LHC, these values can be increased by a factor of order 
3, given the accessibility to partons with smaller momentum fraction in the nuclear wave 
function, or higher plasma temperatures in the hot matter case. We can anticipate that the 
results at the LHC should be qualitatively different since the survival probability decreases 
by about one order of magnitude. 
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We also observe on Fig. [7] that the survival probability decreases as 1/Q 2 S for large values 
of Q s . In the hot matter case, this corresponds to a 1/L behavior in terms of the plasma 



length. This dependence, already derived previously [23J (and in [271, [28( for an SU(2) 
color group), is in contrast with the e~ L decrease usually assumed. Such an exponential 
dependence is valid when the dipole undergoes successful independent scatterings in the 
medium, while the power-law we obtain is due to the coherence of the multiple scatterings, 



and both the elastic and inelastic contributions display this super-penetration feature 



29] 



A similar scaling was also observed when considering the suppressed production of heavy 



quarks off cold nuclear matter [301 ] . 



B. Collisional dissociation of heavy mesons 

While for light hadrons, partonic energy-loss calculations can describe the suppressed 
production of high— px particles in nucleus-nucleus collisions at RHIC, the suppression is 
underestimated in the case of D and B mesons (or rather single non-photonic electrons com- 
ing from the decays) [5J. The discrepancy is due to the small quenching of 5-mesons, which 
dominate the high-p-r single electron yields. However, the standard calculations assume 
that the hard parton hadronizes outside the medium, having fully traversed the region of 
dense nuclear matter, and lost energy via radiative and collisional processes. But for heavy 
mesons, due to the significantly smaller formation times, this assumption does not appear 
to be justified. In fact, it was shown that including collisional dissociation goes in the right 
direction 3]: contrary to calculations that emphasize radiative and collisional heavy quark 
energy loss, collisional dissociation predicts that D and B mesons are suppressed in a similar 
way at transverse momenta as low as pt ~ 10 GeV. 

n 

The framework of Ref. [3[ takes into account the competition between the fragmentation 
of c and b quarks and the medium-induced dissociation of the D and B mesons to evaluate the 
quenching in nucleus-nucleus collisions. The meson survival probability we have computed 
in the previous sections is an important ingredient of their analysis: it determines the 
dissociation time of the meson in the plasma. However the expression used in [3!] ignores the 
color structure of the theory and therefore it would be interesting to redo the analysis with 
our expression, but this is beyond the scope of this work. Note that the time dependence of 
the plasma formation and evolution is easily implemented by varying the parameter Q 2 S . 
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C. Diffractive vector meson production 

An experimental situation where elastic and inelastic processes can be distinguished is 
the diffractive production of vector mesons in deep inelastic scattering: — > VY, where 
A stands for the target nucleus and Y for the final state it has dissociated into. In this 
process, the qq pair that the virtual photon has fluctuated into scatters off the nucleus before 
recombining into a vector meson. While the scattering involves a color-singlet exchange, 
leaving a rapidity gap in the final state, the nucleus can still scatter elastically (Y = A, this 
is called coherent diffraction) or inelastically (i.e. break up, called incoherent diffraction). 
Kinematically, a low invariant mass of the system Y corresponds to a large rapidity gap 
between that system and the vector meson (this also implies that the longitudinal momentum 
of the meson is close to that of the photon, which justifies using the eikonal approximation 
for this process), which makes it possible in principle to keep track of the state of the target, 
and separate coherent and incoherent diffraction. 

The cross-section is peaked at minimum momentum transfer where the elastic scattering 
dominates, but as the transfer of momentum gets larger, the role of the inelastic contri- 
bution increases and eventually it becomes dominant (typically for momenta bigger that 
the inverse nucleus size). The momentum transfer in this process is essentially the trans- 



verse momentum of the vector meson in the final state P\, and as a function of \t\ = P 



± J 



the elastic contribution decreases exponentially while the inelastic contribution decreases 
only as a power law. This important difference was not discussed in this paper where only 
P'x_— integrated quantities are analyzed. It deserves detailed studies which are left for future 
work, such as for instance the numerical analysis of our results and a comparison with data 
from HERA on diffractive vector meson production, with or without proton breakup. The 
case of deep inelastic scattering off a large nucleus should also be studied, and in this case 
the MV model we considered provides a natural framework, and a good starting point to im- 
plement the high-energy QCD evolution. Inclusive and diffractive structure functions have 
been calculated [31]. but vector-meson production has yet to be addressed. At an electron- 
ion collider, when the momentum transfer is small enough for the nucleus to stay intact, 
then it will escape too close to the beam to be detectable; therefore the whole diffractive 
program will rely on our understanding of incoherent diffraction. 
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VI. CONCLUSIONS 

The main technical result of the paper is the 4-point function (|38|) . computed in the MV 
model for cold nuclear matter, but also valid in the GW model for hot matter. This dipole- 
dipole correlator allows to compute the survival probability of a meson propagating in the 
presence QCD matter (TH1) . and the cross-section for the diffractive production of vector 
mesons in deep inelastic scattering off nuclei (fT6l) . In the hot nuclear matter case, only the 
large N c results were explicitely derived, with a diagrammatic approach which gave a more 
physical picture of what are the processes giving the main contribution to the dissociation 
of bound states. It also allowed a clear distinction between elastic and inelastic processes 
with respect to the target. 

The two results for cold and hot matter are given by the same expression when written in 
terms of the saturation momentum of either the nucleus or the plasma (it is more common 
in this case to write Q 2 S = qL) . The medium dependence enters only through this parameter. 
While this allowed us to make general statements, the interplay between the meson wave 
functions and the dipole-dipole correlator should be analyzed numerically. Even though 
experimentally it is not easy to appreciate the difference between elastic and inelastic pro- 
cesses in heavy ion collisions, we consider our result to be interesting from the theoretical 
point of view since we address a realistic situation where the medium can be excited by the 
interaction with the meson. 

We obtained that the survival probability of quarkonia traversing a hot QCD medium 
exhibits the super-penetration feature: it decrease as 1/L for large medium length L. We 
also except our result to be relevant when considering the production of high-pr D and B 
mesons in nucleus-nucleus collisions, as these can be produced within the medium and their 
in-medium dissociation will contribute to the strong suppression observed at RHIC. This 
will be better investigated at the LHC where it should be possible to distinguish signals 
from charm and bottom quarks making it possible to differentiate between the two cases 
in a regime where we can safely rely on the eikonal approximation. Finally, an actual 
comparison of our results for elastic and inelastic processes is more accessible with the 
diffractive production of vector mesons in deep inelastic scattering off nuclei; they correspond 
to coherent and incoherent diffraction respectively which can be derived from our formulae. 
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